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We formulate a theory for entangled imaging, which includes also the case of a large number of 
photons in the two entangled beams. We show that the results for imaging and for the wave-particle 
duality features, which have been demonstrated in the microscopic case, persist in the macroscopic 
domain. We show that the quantum character of the imaging phenomena is guaranteed by the 
simultaneous spatial entanglement in the near and in the far field. 

CN 

A major trend in physics at present is to "push the realm of quantum physics well into the macroscopic world" [0. 
^ In the case of photon systems, quantum effects such as squeezing or twin beams can be found also in the regime of large 
photon number. However, the analysis of fundamental phenomena, such as wave-particle duality, for a macroscopic 
O ■ electromagnetic field requires a spatially multi-mode treatment, as developed in J2|, |3|. 

The field of quantum imaging Q and, especially, the topic of entangled two-photon imaging (EPI) provides an ideal 
framework for such a discussion. The theory of EPI was pioneered by Klyshko j|] who formulated a heuristic approach, 
that stimulated a number of key experiments, especially in the laboratory of Shih Similar experiments in the 
group of Zollinger [Q addressed the discussion of fundamental issues in quantum physics. Recently, the Boston group 
formulated a systematic theory of such phenomena Q . All the papers j^, ^, |?J concerned the regime of single photon- 
pair detection in the parametric downconversion process. In this work we formulate a theory for EPI which holds 
for arbitrary downconversion efficiency, and therefore encompasses also the case in which a large number of photons 
are detected in each pump pulse, as it happens in the experiment described in A key role in our analysis is 
played by the concept of spatial entanglement j^, we introduced previously. Special attention is devoted to a clear 
identification of the features that require the presence of quantum entanglement rather than classical correlations, as 
those of the experiment of Boyd et al. JTTJ . 

We consider a type II x' 2 ' crystal of length l c . We assume an undepleted pump beam, with a Gaussian profile of waist 
w p , and a Gaussian temporal pulse profile of duration r p . In the parametric down-conversion process the signal/idler 
(S/I) photon pairs are emitted over a broad band of temporal frequencies (bandwidth f2o oc l c 1 ) and a broad band 
of spatial frequencies (bandwidth qo oc (A/ c ) -1 / 2 , with A being the central wavelength of the down-converted fields). 
We developed a numerical model, based on the Wigner representation, that simulates the propagation of the three 
waves inside a realistic x*- 2 -* crystal, including the effects of diffraction, spatial and temporal walk-off, and temporal 
dispersion |[2]]. Here we will not describe such a method, but we will focus on some key results for a conceptual 
experimental scheme suited to discuss the wave-particle aspects at a macroscopic level. We also used an analytic 
• *h , approach, valid in the limit of a plane-wave cw pump, where propagation inside the crystal is described in terms of 
[ the unitary transformation ^] : 

of* (q, 0) = Ui (q, Q)af (£ A) + V t (q, H)a? 1 (-q, -O) 

i?3 = S,I (1) 

linking S/I fields at the input with those at the output face of the nonlinear crystal. ai(q,Q) m ^ out are annihilation 
operators of plane-wave modes, q being the transverse wave-vector and O the shift from the carrier frequency. The 
explicit form of the functions Ui and Vi is given for example in || . For brevity, we drop in the following the frequency 
argument from all the formulas (even if we took it into account in our calculations). When the transformation 
equivalent to (Q) in a Schroedinger-like picture is applied to the input vacuum state of the S/I fields, we obtain the 
output entangled state: 

l^) = Il|£ c »®l n '^|n,-g)il , (2) 
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where \n,q)s/i denotes a Fock state with n photons in mode q of the S/I beam, and c n (q) 
{U S (q)Vj(-Q)} n \Us(q)r (2n+1) - Moreover, 

("(g)) 
[1 + M<7)>] 
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where (n(q)) is the average number of photons in mode q. Almost all down-conversion literature is limited to the case 
(n(q)} <C 1, in which state (||) reduces to 

~l[c o (<?)\0,q)s\0,-q)i + \ci(q)\l,Q)s\l,-q)i J] c o(<?i)|0, ?l)s|0, 1 (4) 

In this case, which we refer to as the microscopic case, one detects coincidences of single photon pairs; in application 
to imaging, the image is reconstructed from a statistics over a large number of coincidences. In this paper we focus, 
instead, on the case in which the average photon number per mode is not negligible, so that all the terms of the 
expansion (|^) are relevant (we call it the macroscopic case). In this case the entanglement is with respect to photon 
number, and this model predicts ideally perfect correlations in S/I photon number detected in two symmetric modes 
q and — q JTof ] - 

An interesting analytical limit is that of a short crystal (where diffraction and walk-off of the S/I fields along the 
crystal become negligible), where the output state can be written in the form 

w =n|£ c «(9 = °)i n ' f )si^ a? >4 ; ( 5 ) 

x ln=0 J 

x denotes position in the transverse plane at the crystal exit ( "near- field" ) , and \n, x) is the Fock state with n photons 
at point x. In this limit one has ideally a perfect correlation in the number of S/I photons detected at the same 
near field position. We incidentally note that if only one beam of the two is considered, its reduced density matrix 
is diagonal in the Fock state basis, and corresponds to a thermal statistics with average number of photons given by 

In the more sophisticated numerical model, the finite size of the pump has the effect that, if one idler photon is 
emitted in direction q, its twin photon will travel in the symmetric direction — q, within an uncertainty Sq l/w p , 
which hence represents the uncertainty in the signal transverse momentum when determined from a measurement 
of the idler transverse momentum. On the other hand, due to the finite length of the crystal, twin photons created 
in a single down-conversion process at the same position, are separated by diffraction along the crystal. Hence the 
uncertainty in the position of a signal photon conditioned to the detection of an idler photon at position x is given 
by a coherence length l co h = l/qo ~ (A/ c /27r)2 . Really important for the purpose of imaging is the number of pixels 
that can be resolved in an imaging scheme based on correlation measurements. This number is assessed, both in the 
near and in the far field, by the ratio (w p /l co h) 2 = (qo/Sq) 2 . The simultaneous presence of entanglement both in 
momentum and in position is a fundamental property of the down-converted photons, which, as we will see, plays a 
crucial role in the imaging process. This property persists for a large photon number, a case in which the entanglement 
assures a photon number spatial correlation at the quantum level both in the near and the far field. 
Figure [I] illustrates a compact imaging scheme. The S/I beams are separated by the polarizing beam splitter PBS 
(we assume for simplicity that the distance from the crystal exit to PBS is negligible). In the path of the S beam 
there is an object, which is imaged by a lens on the far-field plane, where it is detected by a single point-like detector 
Ds- An identical lens images the / beam on its detection plane, where it is observed by an array of detectors Dj. 
The distance z between the PBS and the lens, and between the lens and Dj can be varied; we focus on the cases 
z = f and z = 2f, in which we will see that the diffraction pattern (z — /) and the image (z = 2f) of the object can 
be reconstructed by correlation measurements. For definiteness, we discuss the case in which the object is a double 
slit, with a being the width of the two slits and d their distance. For z = /, Fig.|I] a corresponds to the setup of some 
of the experiments in || . 

The object is in the path of the S beam, which is observed by a point-like detector, and no information about it 
can be obtained by direct detection. As a straightforward generalization from the coincidence measurements of the 
microscopic case||, ||, we consider the spatial correlation of the S/I detected intensities. Precisely, we denote with 
Is(xs) and Ii(xi) the intensities detected by Ds and by the array Dj averaged over a detection time td (in typical 
pulsed experiments td w rp), and we introduce the spatial correlation function 

(//(£/)/<?(£<?)) - (Ii(xi))(Is(x s )) + (Sl^SIsixs)) . (6) 

The object information is contained in the correlation of intensity fluctuations G(xi,xs) = (5Ii(xt)5Is(xs)) as a 
function of xj for fixed xs- Toghether with the scheme (a) of Fig.|I], we consider the alternative scheme (b) in which, 
conversely, the S beam is detected by an array, and / by a point-like detector. Such a scheme was analysed in [jjj in 
the microscopic case. In the macroscopic case the image is provided by G(xi, xs) as a function of xs for fixed a?/. 
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FIG. 1: (a)Imaging scheme. P = pump beam, x' 2 ' = type II crystal, S = signal field, I = idler field, Ds,Di = detectors, L = 
lens with focal length /, NF = near field, FF = far field, PBS = polarizing beam splitter, (b) scheme for the discussion of 
fundamental aspects. 



Let us first consider the case z = /, in scheme (b). If the / field is not detected, there is no possibility of observing 
the interference fringes by direct measurement of field S alone, unless the object is contained in a coherence area, i.e. 
d <C l C oh- In the microscopic case it was argued [[l] that in principle one could detect the I photon, and obtain "which- 
path" information on the S photon, and this is enough to cancel the fringes. More in general, we argue that, since the 
S beam alone is in a incoherent thermal mixture, the interference fringes are not visible due to the lack of coherence. 
However, in order to make fringes visible, it is enough to condition the S beam measurement to a measurement of the 
/ beam by a single point-like detector. In the microscopic case the fringes are observed via coincidence measurements, 
as explained in fl , because detection of the / photon in the far field determines the S photon momentum before the 
double slit, due to momentum entanglement, providing a quantum erasure of any which-path information. In 
the general case, the basic mechanism is the S/I far field intensity correlation, and calculations performed with the 
analytical model (0) show that 



G(xs,xi) oc 



T q 



2tt 

a/ 



(x s + Si) 



(7) 



where T(q) is the Fourier transform of the transmission function T(x) describing the object. Under the conditions 
a }t Icoh, d < w p , the entire interference-diffraction pattern is visible with good resolution. The result (pi) is symmetric 
with respect to xs and xj, hence the same pattern appears in the imaging scheme (a). Figj^b shows the result of a 
1-D numerical simulation of the pattern reconstruction via intensity correlation function. A statistical average over 
a reasonable number of pump shots was enough, because t p was on the same order of magnitude as the amplifier 
coherence time r co ^ = 1/CIq. Our calculations show that when rrj S> t co /j, the first term at rhs of Eq. (^) becomes 
much larger than the second term, that contains all the information about the object, at the expenses of the visibility. 
Consider now scheme (b) in the z = 2f case, in which the Di detector lies in the image plane with respect to the 
object, and the measurement exploits the S/I spatial correlation in the near-field. In the microscopic case fringes 
are not visible because the detection of the / photon in the near field, due to position entanglement, provides perfect 
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FIG. 2: Numerical simulation of the experiment in Fig. 2a with z — f. Parameters are those of a 4 mm BBO crystal, with 
w p — 332pm, t p — 1.5ps (r co h = 0.87ps, l co u = 16.6 (im), a = 17pm, d = 104pm. (a)Mean intensity of the S/I beams after 
10000 pulses; (b) Solid line: correlation G(xi,xs) as a function of xi after 10000 pulses; dashed line: plane-wave result of 
Eq.(§. x = Xfq /(2n) 
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FIG. 3: Numerical simulation of the experiment in Fig. 2a with z = 2/. Parameters as in figure ^ The figure shows G(xi,x,s) 
as a function of xi after 10000 pulses 



"which path" information about the S photon jjj. Our general result is that, again, for a > l co h : d < w p 



Us (^) Vl (.«-) 



(8) 



In scheme (b), where Xi is fixed, there is no information about the object. However, in scheme (a) (xs fixed)the 
object image can be reconstructed via the correlation measurement. Hence, by only changing the optical setup in the 
path of the idler, which does not go through the object, one is able to pass from the diffraction pattern to the image 
of an object. This result is confirmed by our numerical simulation shown in Figj^. 

In order to assess the quantum nature of the phenomena observed in the imaging scheme (a), the key question is 
whether these results can be reproduced by using a "classical" mixture, instead of the pure entangled state (|]|]). It 
is natural to focus on the two mixtures 



W = n g -{E,T=o \cn{q)\ 2 \n,q)s\n,-q)n(n,q\ s (n,-q\} 



(9) 
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W = Us {En= K(0)\ 2 \n^s\n^H(n,x\s(n,x\} (10) 

Mixture (§) preserves the local S/I spatial intensity correlation in the far field, while the intensity correlation function 
is completely delocalized in the near field. By following the same notation of |?J , we indicate by hj (xj , x) the linear 
Kernel describing propagation throught the imaging setup of beam j = I, S; we introduce their Fourier transforms 
hj (xj , q) describing how a q component of the j beam at the crystal exit face is transformed into the field at point Xj 
at the detection plane. With mixture (0) we obtain 



G(x s ,x I )= d^hsiSs^h^-mUsimi-^ 1 . (11) 
The optical setup in the S beam arm is fixed, and 

hs(x S ,q) = -^-J C^-xs - q) (12) 



In the z = f configuration of Fig. 2a, \hi(xi, — q)\ 2 oc 5(jjXi + q) and we obtain the same result of Eq.(@). Hence 
fringes are visibile with the classical mixture (||) in the same way as with the pure EPR state @. However, for 
z = 2/, hj(xi,—q) oc exp (ixi ■ q), and the correlation function is constant with xi; thus in this case the scheme gives 
no information at all about the object. Conversely, the mixture (|l0| ) preserves the S/I local intensity correlation only 
in the near field. Not surprisingly, in this case the z = 2/ scheme (a) provides the image of the object, as with the 
pure state, but in the z — f case the fringes are not visibile. The key point is that only the pure EPR state (||J|) 
displays S/I spatial correlation both in the near and in the far field. This analysis is not in contrast with the basic 
conclusion of Ref. |ll| , that the result of each single experiment in EPI can be reproduced by a classically correlated 
source. Here we argue that only in the presence of quantum entanglement the whole set of results illustrated in Fig. 2b 
and 3 can be obtained by using a single source, and by keeping the optical setup in the signal beam arm fixed. 
In conclusion, we formulated a theory that encompasses both the microscopic (single photon pair detection) and the 
macroscopic (multi-photon detection) case. Our results show that the imaging and wave-particle duality phenomena, 
observed in the microscopic case, persist in the macroscopic domain, and indicate a possible experiment that is able 
to discriminate between the presence of quantum entanglement or classical correlation in the two beams. Clearly, 
there is a pratical limit in the macroscopic level that can be attained preserving such phenomena. In order to increase 
the number of down-converted photons, usually either the pump beam is more focused (w p is decreased), or the 
crystal length l c is increased. However, when the condition w p — l co h oc \J Xl c is reached, the resolution of the spatial 
entanglement in the near and far field, as well as in the entangled imaging, is completely lost. 
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